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due to their attractive features of being renewable in nature and reducing the net CO2 emissions. Biofuels 
have been used in conventional diesel and gasoline engines either as neat fuels or as supplements. 


Key Mon Fortunately, a relatively new combustion concept for internal combustion engines, namely homo- 
Biofuels z SE ! g à 

HCCI geneous charge compression ignition (HCCI) combustion, has been evolved in parallel to the biofuel 
Simulation research. HCCI combustion seems to be able to take advantage of the diverse properties of biofuels, since 
Ethanol in this combustion mode ignition is not externally instigated, but relies on the compression and sub- 
Hydrogen sequent autoignition of a fuel-air mixture. This fact allows the utilization of different fuels or blends 
DME thereof, in order to regulate the ignition point and provide adequate operation under diverse operating 
DEE conditions. 

Biodiesel This study provides an overview of existing simulation models for the simulation of biofueled HCCI 
Methyl butanoate combustion. Simulation models aid and supplement the experimental research conducted on HCCI com- 


Methyl decangate bustion, providing a fundamental insight into the physicochemical parameters affecting performance and 


emissions formation. The simulation models include single-zone models, multi-zone models, probability 
based models, and multi-dimensional models in order of complexity. The vast majority of these models 
implement chemical kinetics to simulate the combustion process, not only due to the inherent depen- 
dence of HCCI combustion on the physicochemical properties of the fuel, but also due to the sometimes 
complex chemical structure of the biofuels, which include esters, ethers and alcohols. The reaction paths 
for these homologous series are quite different from the conventional hydrocarbons used to simulate 
conventional fuels, and provide the ground for current and future research work. 

© 2011 Elsevier Ltd. All rights reserved. 
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Nomenclature 

B cylinder bore (m) 

h molar specific enthalpy (J/kmol) 
i any zone 

k species index 

K number of species 
n moles 

P pressure (N/m?) 

S cylinder height (m) 
T temperature (K) 

t thickness (m) 

U internal energy (J) 
V volume (m?) 

Z number of zones 


Greek symbols 


À air-fuel equivalence ratio (lambda) 
(o fuel-air equivalence ratio (phi) 
Subscripts 

chem chemical 

cyl cylinder 

HT heat transfer 

i zone index 

in inlet or initial 

w wall 

Superscripts 

in inflow 

out outflow 

tr transferred 


Abbreviations 


a after 

b before 

BDC bottom dead center 

CA crank angle 

CAD crank angle degrees 

CO carbon monoxide 

CFD computational fluid dynamics 
CI compression ignition 

CR compression ratio 

DI direct injection 


DTBP di-tertiary butyl peroxide 
EDC eddy dissipation concept 


EGR exhaust gas recirculation 
EtOH ethanol 

EU European union 

EVC exhaust valve closing 

EVO exhaust valve opening 

FID flame ionization detector 
HC (unburned) hydrocarbon(s) 


HCCI homogeneous charge compression ignition 


HRR heat release rate (J/deg) 


HTR high temperature reaction 
HTHR high temperature heat release 
IA injection angle 


imep indicated mean effective pressure (bar) 
ISFC indicated specific fuel consumption (g/kWh) 


IT injection timing 

IVC inlet valve closing 
IVO inlet valve opening 
LES large eddy simulation 


LTR low temperature reactions 
LTHR low temperature heat release 
MB methyl butanoate 

MD methyl decanoate 


MeOH methanol 

MON motor octane number 
MRG methanol reformed gas 
MSW municipal solid waste 


MZ multi-zone 
NOx nitrogen oxides 
NVO negative valve overlap 


OHC oxygenated (unburned) hydrocarbons 
PDF probability density function 


PI port injection 

PRF primary reference fuel 
RES renewable energy sources 
RON research octane number 
SRM stochastic reactor model 
SI spark ignition 

SZ single-zone 

TC turbocharged 

TDC top dead center 

toe tonnes oil equivalent 


UHC unburned hydrocarbons 


1. Introduction 


During the last decades two main concerns have arisen in the 
energy generation and transportation sectors, namely the increas- 
ing energy demand coupled with the limited availability of oil 
reserves and the detrimental environmental effect produced by 
their exploitation to deliver energy and motion. Nuclear energy 
could be used to gradually displace fossil fuel based energy pro- 
duction, but the late accident in the Fukushima nuclear power 
plant has raised many concerns related to the hazards of nuclear 
power plants; it seems that the efforts toward the establishment 
and the public acceptance of nuclear energy have yielded some- 
what to skepticism. In an effort to gradually displace fossil fuels, 
gross inland energy consumption! of renewable energy sources 
(RES) has increased in the European Union (27 countries) in the 
decade 1999-2009 from over 5% to approximately 9% [1]. How- 
ever, this increase in the share of renewable energy consumption 
seems to be mainly due to the reduction in the solid fuel energy 
consumption share, since the nuclear energy consumption share 
(~13-14%) and the combined share of petroleum and natural gas 
(~61%) have remained fairly constant (Fig. 1). 

The use of fossil fuels has increased the energy dependence of 
many countries around the world, since oil production is extremely 
localized with the Middle East possessing over 55% of the proved 
oil reserves world wide as of 2009 [2,3]. This impacts all countries 
and citizens around the world, since the oil prices are reflected 
on virtually every commodity. Moreover, it creates and sustains 
energy insecurity among countries with high energy dependence. 

Another detrimental effect of the use of fossil fuels is related to 
the environmental impact from their use for energy production and 
transportation. The fuel combustion in internal combustion engines 
produces - among others - CO2, a major greenhouse gas that is 


1 Gross inland energy consumption refers to the total energy demand of a country 
or region. It covers energy consumption by the energy sector itself, distribution 
and transformation losses, final energy consumption by end-users and statistical 
differences. It does not include energy provided to maritime bunkers [1]. 


1590 N.P. Komninos, C.D. Rakopoulos / Renewable and Sustainable Energy Reviews 16 (2012) 1588-1610 


100% 


Renewable Energy Sources 


90% 
80% 
70% 
60% 
50% 
40% 
30% 
20% Total Petroleum Products 


10% 


0% 


1999 2001 2003 2005 2007 2009 
Year 


Fig. 1. Gross inland energy consumption share by fuel in the EU (27 countries). Gas: 
almost exclusively natural gas. Renewable energy sources: biomass and renewable 
wastes, hydro, geothermal, wind, solar, wave, tide and ocean energy. 


Data source: Eurostat (online data codes: nrg_100a, nrg101a, nrg_102a, nrg_103a, 
nrg.104a, nrg.1071a) [1]. 


known to contribute to the global warming effect. The main issue 
is not the production of CO; itself, since the combustion of most 
fuels leads to its production, but the fact that in the case of fossil 
fuels the CO; emitted cannot be reconverted to the original fuel by 
any short-term, natural - and therefore energy neutral - process. 

In view of the first issues analyzed, i.e. the energy dependence 
of most of the countries, the instability of the oil prices and the 
increasing demand of fossil fuel energy, alternative fuels have been 
sought after. Among those, biofuels have received increasing atten- 
tion due to their additional attribute of being renewable. The term 
"biofuel" refers to any fuel that can be produced from biomass, 
such as starch, sugars, lignocellulosic material, animal fat, oleagi- 
nous plants, etc. The main advantage of a biofuel is that it can be 
reconstructed by natural processes using the CO; emitted from its 
combustion. This is part ofthe photosynthetic process ofthe plants, 
which use sunlight and CO; to build their cells. Thus, biofuels have 
the potential to be CO; neutral. This is depicted in Fig. 2, in which 
a comparison is made between the biofuel energy system and the 
fossil fuel energy system. 

Initially, feedstock such as barley, wheat, corn, sugar beat, sugar- 
cane, etc. was - and still is - used for ethanol production; oleaginous 
plants such as soybeans, castor beans, rapeseed, Jatropha and palm 
were used for biodiesel production. However, the feedstock utilized 
for the production of these "first generation" biofuels originates 
from edible crops. This led to an unfavorable effect related to food 
security in terms of both land use and food prices. With over 1/8th 
of the population worldwide undernourished as estimated in 2010 
[5,6], using part ofthe edible crops to produce biofuels did not seem 
a socially just choice. Added to the land use was the issue of rising 
food prices due to the increased demand, which was instigated by 
the alternative use of food commodities for fuel production [2]. 

The food security issue raised concerns related to the choice of 
biofuel production routes and shifted research toward biofuel pro- 
duction from non-edible feedstock, mainly lignocellulosic material 
or "plant biomass". The term "second generation" biofuels is used 
to depict this shift in the feedstock used, but not the final fuels 
produced. The lignocellulosic material is abundant in plants and 
mainly consists of polysaccharides [7]; biofuels can be obtained 
from this feedstock by hydrolysis and fermentation, for ethanol 


production, and by gasification to produce liquid hydrocarbons, 
biodiesel, bioDME and bio Synthetic Gas. Additionally, algae seem a 
promising feedstock for the production of second generation biofu- 
els and especially diesel fuel. Their advantages are their high growth 
rates, tolerance to diverse environmental conditions, high CO; trap- 
ping and fixation ability, high oil content - sometimes exceeding 
80% of their dry weight - and the high microalgae yield per hectare, 
which is almost an order of magnitude higher than the palm oil 
yield [7]. Owing to the fact that other non-fuel products can also be 
obtained from biomass, an integrated approach is being researched 
to fully exploit the biomass available and produce all possible prod- 
ucts from biomass, including biofuels and a platform of chemicals, 
in biorefinery systems [4,7]. The production of second generation 
biofuels is not yet cost effective since there are a lot of technical 
barriers to be overcome [7]. 

In view of the attractive attributes of biofuels, efforts have been 
made to introduce them into the energy consumption systems. In 
Fig. 3 the final energy consumption? of RES in EU is shown for the 
year 2009. It can be seen that the main sector in which biofuels — 
especially liquid ones - are currently used, is the transport sector 
and particularly road transports. Fig. 3 also indicates that liquid bio- 
fuels - including bioethanol, biomethanol, biodiesel, bioDME and 
other - are currently the gateway for the penetration of RES into 
the transport sector. The share of liquid biofuels in the transport 
sector and especially road transports has increased during the last 
decade as shown in Fig. 4. 

These data highlight the importance of biofuels and justify the 
research effort toward their utilization in conventional (CI and SI) 
engines as well as in Homogeneous Charge Compression Ignition 
(HCCI) engines. In Diesel engines fuels such as n-butanol, ethanol, 
biodiesel, DME, either neat or in blends, have been tested and/or 
simulated [8-16]. Some of these biofuels, such as esters and DME, 
have an inherent tendency to reduce soot formation due to their 
chemical structure; the tendency to form soot precursors decreases 
as the oxygen content in the fuel molecule increases and/or the 
number of C-C bonds decreases [16]. This is due to the fact that 
the oxygen bound carbon atoms tend to form CO or CO; rather 
than participate in soot formation reactions [17]. In spark ignition 
engines hydrogen, biogas, ethanol, DME have been tested and/or 
simulated [18-21]. 

Apart from the conventional combustion modes a relatively new 
concept has been developing during the last two decades, i.e. homo- 
geneous charge compression ignition combustion (HCCI). In HCCI 
combustion a mixture of fuel and air is introduced into the combus- 
tion chamber early during compression or during the intake stroke. 
The mixture is compressed and autoignites at the point where the 
local thermo- and fluid-dynamic conditions are favorable [22]. The 
main advantage of the HCCI concept relative to conventional com- 
bustion modes is the significant reduction of NO, emissions and 
the absence of soot in the exhaust. Moreover, the fast combustion 
usually encountered and the unthrottled operation are prone to 
yield high indicated thermal efficiency. However, HC and CO emis- 
sions are usually high due to bulk quenching at low loads, and wall 
quenching and crevice flows at higher loads [23]. 

Added to these unfavorable attributes of HCCI is the inability 
to directly control the combustion phasing in HCCI engines. This is 
a major obstacle preventing the implementation of the HCCI con- 
cept. HCCI ignition relies heavily on chemical kinetics and largely 
depends on the fuel used. In an effort to control ignition and the 
combustion rate, thereby expanding the HCCI operation over a wide 
range of load-engine speed, various means have been used, such 
as variable compression ratio, variable valve timing, variable inlet 


? Final energy consumption is the total energy consumed by the end-users, such 
as industry, transports, households, services, etc. [1]. 
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Bioenergy system Fossil reference energy system 


Fig. 2. Comparison of bioenergy and fossil fuel energy systems [4]. 


temperature, internal or external EGR, direct fuel injection, etc. [23]. conditions [23-25]. The main concept is that fuels with different 
A promising method for the combustion control seems to be the autoignition tendencies can be blended at varying proportions to 
use of alternative fuels and fuel blends, depending on the charac- regulate the ignition point at various load-speed regions [26]. This 


teristics of the engine (compression ratio, etc.) and the operating is a fortunate evolution since the fuels used in the blends may 
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Data source: Eurostat (online data codes: nrg.1071a, nrg.1073a) [1]. Data source: Eurostat (online data codes: nrg. 100a, nrg.1073a) [1]. 
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as well be biofuels. Thus, hydrogen [27], ethanol [28-32], ethers 
[33-37] and biodiesel [38] have been used in HCCI engines, either 
as neat fuels or in blends. High cetane number fuels (such as DME) 
can be used as ignition promoters in the blend and high octane 
number fuels (e.g. ethanol, methanol, etc.) can be used as ignition 
suppressors. With the use of closed or open loop control systems, 
the resulting fuel blend can provide the means to control HCCI igni- 
tion over a wider range than the one provided by the individual 
neat fuels [25]. This has already been demonstrated in HCCI engines 
using blends of n-heptane and isooctane [26]. 

As is the case with conventional combustion types (SI and CI), 
the HCCI combustion mode can be studied using simulation mod- 
els of various degrees of complexity. Simulation models describe 
the fundamental processes occurring within the combustion cham- 
ber and - if accurate enough - can help toward understanding 
the effects of individual operating parameters on HCCI combus- 
tion. Such an analysis is not always feasible in experiments, since 
it is usually difficult to isolate the effect of a single parameter 
- be it fuel properties, inlet conditions, etc. - without directly 
or indirectly affecting another parameter. Simulation models can 
provide insight into the fundamental processes occurring during 
combustion and/or during the stage of gas exchange processes, 
depending on the sophistication level of the model. Features such 
as combustion chemistry, heat and mass transfer can be studied 
so that the crucial factors affecting them can be assessed. Espe- 
cially in the case of biofuels, which have not been numerically 
studied extensively until recently, there is a lot to be explored 
as regards combustion chemistry, emissions formation, trans- 
port properties, etc. Combustion chemistry affects significantly 
the mixture autoignition, combustion rate and emissions forma- 
tion. Transport properties affect heat and mass transfer within the 
combustion chamber, which determine the temperature field or 
temperature-mass distribution within the combustion chamber. 
Moreover, in the case of direct injection HCCI, in which the fuel is 
directly injected into the combustion chamber, the fuel transport 
properties must be adequately modeled to simulate the forma- 
tion, break-up, vaporization and mixing of the fuel jet with the 
air. Such phenomena determine the mixture stratification and the 
possibility of fuel jet impingement on the combustion chamber 
walls. 

HCCI simulation models are divided into single-zone, multi- 
zone and computational fluid dynamics (CFD) models, in order of 
increasing complexity. In single-zone models the main assumption 
is the uniformity of the charge properties within the combus- 
tion chamber. Multi-zone models consider multiple regions within 
the combustion chamber, each of which is considered an open 
thermodynamic system of uniform properties. CFD models are 
the most sophisticated ones, providing detailed description of 
the in-cylinder gas motion and thermodynamic state on a local 
level. However, this increased level of sophistication is not with- 
out compensation, in terms of increased computational time and 
effort. 

In the following sections single-zone, multi-zone and CFD mod- 
els, which have been constructed for biofueled HCCI combustion 
are presented. Since the vast majority of HCCI simulation models 
use chemical mechanisms of various levels of detail to describe 
the chemical kinetics of the ignition and combustion processes, a 
separate section is devoted to the use of chemical kinetics. The 
HCCI simulation models presented in the following sections are 
included in this study regardless of whether they produced specific 
biofueled-oriented results. The selection criteria were considered 
satisfied by the mere fact that biofueled HCCI combustion was 
simulated. Upon consideration of the relatively short simulation 
history of the HCCI engine, simulation data for biofueled HCCI 
engines have their own obvious merits and provide a reference for 
future studies nonetheless. 


2. Chemical kinetics modeling - reaction pathways 


Almost all of the simulation models used in HCCI combus- 
tion modeling and all simulation models presented herein include 
chemical kinetics. This is the result of the decisive role chemical 
kinetics play in HCCI combustion. Since no external means are uti- 
lized to instigate ignition, the latter relies on the reaction kinetics of 
the fuel. Ignition commences where the thermodynamic conditions 
of pressure, temperature and species concentration are favorable. 
Moreover, the specific characteristics of the reaction paths and 
intermediate species also affect the formation of pollutants such as 
NOx, CO and unburned HC. These issues provide the momentum for 
the ongoing research oriented toward building chemical reaction 
paths, which are able to describe in various levels of detail the oxi- 
dation characteristics of fuel-air mixtures under various operating 
conditions. Compared to the diversity of the biofuels used, which 
include hydrogen, alcohols, ethers and esters, the chemical mech- 
anism research has a short history. The difficulties encountered 
involve the identification of the primary and secondary reaction 
pathways as well as the reduction of the resulting chemical mech- 
anisms. This reduction is paramount if the chemical mechanisms 
are to be used in simulation models, since the simulation time 
increases non-linearly with the increase of species and reactions. 
Such considerations are especially important when using the chem- 
ical mechanisms in CFD models, in which the computational time 
cost is already aggravated by the solution of the Navier-Stokes 
equations. 

As regards the chemical mechanisms used in modeling HCCI 
combustion, care should be exercised during the choice of the spe- 
cific mechanism, since the various existing mechanisms do not 
necessarily provide similar results even for well-known fuels such 
as isooctane [39]. Therefore, comparisons of ignition timings and 
global combustion rates must be interpreted with caution, espe- 
cially when differences between fuels are assessed [40]. This also 
holds true for the comparison of NO, emissions, which are known to 
be sensitive to peak temperature variations and consequently to the 
ignition timing predictions. Thus, it is essential for the researcher 
to investigate the conditions under which the mechanisms were 
developed and the experiments used for their validation (shock 
tube, flow reactor, jet stirred reactor, rapid compression machine, 
laminar flame, HCCI engine, etc.). Preferably, the mechanisms have 
to be validated at low equivalence ratio and high pressure since 
these conditions are relevant to HCCI combustion in engines [39]. 
Furthermore, it has to be kept in mind that the fuel oxidation chem- 
ical mechanisms usually refer to surrogates of actual fuels, which 
cannot describe in exhaustive detail all the actual fuel properties. It 
is obvious that even the detailed mechanisms have their limitations 
and their validation cannot be exhaustive, considering the diverse 
operating conditions encountered in HCCI engines (high EGR dilu- 
tion, fuel injection in hot residuals during negative valve overlap 
HCCI, etc.). 

Table 1 includes chemical reaction mechanisms for various bio- 
fuels found in the literature. Some of these mechanisms have not 
been used in the simulation studies presented herein, but are avail- 
able for future research. However, the vast majority of the oxidation 
mechanisms included in Table 1 has been implemented in sim- 
ulation models, which are described in detail in the subsequent 
sections. The citations describing each mechanism have been num- 
bered consecutively, while the citations of the papers in which 
these mechanisms were used have been numbered in order of 
appearance in the text. 


3. Single-zone modeling 


Single-zone models are often used in HCCI modeling due to 
the short computational time required for the simulation. In the 


Table 1 


Available reaction mechanisms for biofuel combustion simulation. 


# Number of Source Ref. Used in Ref. NO,? Fuel Simul. model Engine type’ CR Pint (bar) Tint (K) Load ($)* Engine Validation! Comments 
species/number speed 
of reactions (rpm) 

1-Butanol (BuOH) 

1. 181/1703 and Dagaut et al. - w n-heptane 
573[2701 [41] 

2. 234/1399 Black et al. [42] - 

3. 251/1990 Dagaut and - w gasoline sur. 

Togbé [43] 
4. 1046/4398 Saisirirat et al. Saisirirat et al. - BuOH/n-heptane SZ 4S 16 - 353 0.3 1500 QLT 63% n-heptane 
[44] [100] 37% BuOH 
(vol.) 

Diethyl ether (DEE) 

5. 112/484 (+DME, Macketal.[45] Mack et al. [45] DEE/EtOH/DTBP SZ 4S 1625 17 380-455 0.3 1800 QLT 
EtOH) 

6. 215/1051 Yasunagaetal | - 

[46] 

Dimethyl ether (DME) 

7. 23/23 Yamada [47] - 

8. 26/28 Yao et al. [48], Huang et al. - DME CFD 4S 17.0 1.1 375 0.01-0.14 1400 QNT 

Liang et al. [49] [138] 

9. 27/35 (+MeOH) Liang [50] Chen et al. [82] DME-MeOH CFD 4S 17.5 1.02 375-380 0.224, 0.226, 1400 QNT DME PI, MeOH 
30/44 [133] 0.605 DI. 
(+MeOH + NOx) 

Yaoetal[137] - DME/MeOH CFD 4S 17.0 1.16 375.6 0.238-0.5 1400 QNT 

10. 28/45 Kim et al. [51] - 

11. 78/336 Curran et al. Mack et al. [45] See DEE 

[52] 
Yao et al. [48] - DME SZ 4S 17 1 348 0.215, 0.095 1400 QNT 
Mack et al. - DEE/EtOH SZ 4S 16.25 1.7 - 0.3, 0.4 1800 QLT 
[101] 
Yaoetal[103] - DME/MeOH SZ 4S 17 1 350 dome - 0.17 1400 QNT 

meon — 
0.063 
Yaoetal[104]  - DME/CH4 SZ 4S 17 1 348 0.297 1400 No 
Ogawa et al. - DME/MeOH-H20- SZ 4S 18.8 1 293 0.16, 0.25 1200 No 
[105] EtOH-H»; -CH4 
Yamada et al. - DME/MeOH-05 SZ 4S 8 1 378-388  0.12-0.42 600 QNT 
[108] (bme) 
12. 79/351 Fischer et al. Hoffman and [83,84] DME SZ 4S 16 1.07-1.49 300 0.4 2300 No 

[53], Curran Abraham [97] 

et al. [54], 

Kaiser et al. 

[55] 
Shudo and - DME/H2 SZ 4S 9.7 - Ambient 0.27-0.35 1000 QLT 
Yamada [106], 
Shudo et al. 
[107] 
Kim etal. [120] [84] DME CFD 4S 17.8 320 8 mg/stroke 1500 QNT IT 80°-0° bTDC 
Park [139] [84] DME CFD 4S 17.8 1-3 300-460  0.2-2.0 1500 No 
Kong [142] [84] DME/MeOH CFD 4S 17.7 1.1 340 0.175-0.435 960 QNT 

13. 97/457 Chun-lan et al. Chun-lan et al. [85] DME CFD 4S - 1 350 0.214 1500 QNT Emissions 
(Detailed) 36/73 [56] [56] validation 
(Reduced) 
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Table 1 (Continued) 


# Number of Source Ref. Used in Ref.? NO,P Fuel Simul. model Engine type‘ CR Py (bar) — Tin’ (K) Load ($)* Engine Validation’ Comments 
species/number speed 
of reactions (rpm) 
Ethanol (EtOH) 
14. 38/228 Röhl and Peters Komninosand [86] EtOH MZ 4S 21 1 ~393 0.1-0.28 994 QNT 
[57] Rakopoulos 
[40] 
Komninosand [86] EtOH MZ 4S 21 1 ~393 0.15-0.26 994 QNT 
Rakopoulos 
[113] 
15. 39/238 Li et al. [58] Sjoberg and - - - - - - - - - - Examined, not 
Dec [118] used 
16. 43/235 Saxena and Viggiano and - EtOH CFD 4S 21 1 363-393 0.16-0.26 1000 QNT 
Williams [59] Magi [145] 
17. 57/383 Marinov [60] Szybist [95] - EtOH SZ 4S 7-29 1 325,500-550 1, ~0.25 900, 2000 No NVO 
Mack et al. [45] See DEE 
Ng and [84] EtOH SZ 4S 15 1 380-660 0.2-0.5 1000 No EGR, neat EtOH 
Thomson [92] or Reformed 
Martinez-Frias [84] EtOH (wet) SZ 4S 16 3 420 (ivc) 1 (deton) 1800 No TC 
et al. [93] 
Yang et al. [96] - EtOH SZ 4S 13 - - - 1400 No NVO 
Mack et al. - EtOH/DEE SZ 4S 16.25 1.7 - 0.3, 0.4 1800 QLT 
[101] 
Ogawa et al. See DME 
[105] 
Komninos and [86] EtOH MZ 4S 21 1 ~393 0.15-0.26 994 QNT 
Rakopoulos 
[113] 
Sjoberg and - - - - - - - - - - Examined, not 
Dec [118] used 
Yu et al. [143] - EtOH CFD 4S 17 1 354, 400 0.303 1200 QNT 
Joelsson et al. - EtOH CFD 4S 17.247 0.9-1 354, 388 0.303 1200 QNT 
[144] 
18. 58/310 Curran et al. Sjoberg and - EtOH MZ 4S 12.8 0.84-1.81 418-443 0.25-0.6 500-2400 QNT 
[61] Dec [117] 
Sjoberg and - EtOH MZ 4S 12.8 1 427 0.4 1200 QNT 
Dec [118] 
19. 112/484 Curran in Mosbach etal.  - EtOH/DEE SRM -PDF 4S 16.25 1.18 400-420 0.3, 0.4 1800 QNT 
Mosbach et al. [62] 
[62] 
20. 167/1591, Dagaut and - w n-heptane 
564/2589 Togbé [63] 
21. 235/1866 Dagaut and - w gasoline surr. 
Togbé [64] 
22. 1046/4398 Saisirirat et al. Saisirirat et al. - EtOH/n-heptane SZ 4S 16 - 353 0.3 1500 QLT 63% n-heptane 
44] [100] 37% EtOH (vol.) 
Ethyl acetate (Ethyl Acet.) 
23. 23/142 Gasnot et al. Cotino and - Ethyl Acet. CFD 4S 13.3 1.35 boost 462 1.5-4 bar 1200 QNT Validation with 
65] Jeanmart [150] imepg iso-octane 
Hydrogen 
24. 10/19 Mueller et al. - 
66] 
25. 10/19 Marinov et al. Goldsborough - H2 SZ 2S 3-38 1.5 300 0.15-0.6 - No Free piston 
67] and Van engine 
Blarigan [89] 
Komninosetal. [86] H2 MZ 4S 18-20 1-2 390-410 0.2-0.4 1000-3000 No 


[112] 
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Table 1 (Continued) 


# — Number of Source Ref. Used in Ref? NO, Fuel Simul. model Enginetype^ CR Pin (bar) Tin? (K) Load ($)* Engine Validation’ Comments 
species/number speed 
of reactions (rpm) 

26. 11/19 O'Connaire - 

et al. [68] 
27. 11/23 Westbrook and Noda and [83] H2 MZ 4S 20 345-410 0.1-1.0 1600 No 
Dryer [69] Foster [111] 
28. 11/25 Kee et al. [70] Fiveland and - H2 SZ AS 15 1.2.1.5 400-800 0.15-0.4 1500 No 
Assanis [88] 
29. 12/30 Li and Karim Liu and Karim [71,72] H2 CFD 4S 16-20 0.85 280-320,515 05 800, 900 QLT 
[71], Liu [72] [149] 

Methanol (MeOH) 

30. 20/84 Li et al. [73] - 

31. 22/89 Held and Dryer Yao et al. [103] See DME 

[74] 
32. 26/84 Westbrook Ogawa et al. See DME 
et al. [75,76] [105] 
33. 52/326 Lindstedt and Yang et al. [96]  - MeOH SZ 4S 13 - - - 1400 No NVO, 32 
Meyer [77] Species/167 
Reactions 

Methyl butanoate (MB) 

34. 41/150 (neat) Brakora et al. Brakora and [87] MB + n-heptane SZ 4S 18 1.1 300-400 0.298-1.0 1200 No 56 Species/169 
53/156 (w [78] Reitz [98] Reactions 
n-heptane) 

Um and Park [84] MB + n-heptane CFD 4S 17.8 1-3 boost 300-460 0.2-2.0 1500 No 
[152] 
35. 275/1545 Dooley et al. - 
[79] 
36. 295/1498 Gail et al. [80] Szybist et al. [84] MB * n-heptane SZ 4S 10.5 1 433-498 0.38-0.48 1800 QNT Tin variable 
[38] 
Methyl decanoate (MD) 
37. 3036/8555 Herbinet et al. Hofmann and [83,84] MD SZ 4S 16 1.16-1.65 310,315 0.4 1670 No 
[81] Abraham [97] 


QLT, qualitative validation-mostly validation of predicted trends; QNT, quantitative validation-includes direct comparison of simulated quantities to experimental measurements; NUM, numerical validation-comparison of the 
HCCI simulation results using a reduced mechanism against predictions from a detailed mechanism. 

? The source mechanisms may have been modified for their use in each study. Moreover, some of the mechanisms may have been updated since the publication of each study. Thus, the number of species and reactions are 
subject to change. Symbol "-" used in this column indicates mechanisms, which are not used in the studies presented herein. 

b Symbol "-" in NOx: information not available in the study. 

* AS: four stroke engine, 2S: two stroke engine 

4 Subscript "in" refers to inlet or initial conditions 

* Load is expressed in terms of the fuel-air equivalence ratio, unless explicitly stated otherwise. 

f Validation refers to the corresponding HCCI simulation model and not to the reaction mechanism(s). 
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Fig. 5. Schematic of the First Law of Thermodynamics for single-zone models. 


single-zone (or zero-dimensional) model the engine cycle is mod- 
eled using the First Law of Thermodynamics, taking into account 
the work absorbed/produced by the piston and - in some cases - 
the heat exchanged between the charge and the combustion cham- 
ber walls, as shown in Fig. 5. Mass blow-by is usually neglected as 
shown in the figure. 

The thermodynamic (pressure, temperature, composition, etc.) 
and transport properties (viscosity, conductivity, etc.) of the mix- 
ture are considered uniform. Since any temperature stratification 
within the combustion chamber is neglected, the single-zone 
models tend to overestimate the global heat release rate and, con- 
sequently, the pressure rise rate and the maximum combustion 
pressure. Moreover, since the mixture composition is also consid- 
ered uniform, the single-zone model is incapable of describing the 
HC and CO formation processes, unless these are due to global par- 
tial combustion of the mixture; under normal HCCI combustion 
conditions, the formation of these pollutants is believed to occur in 
the relatively low temperature regions of the combustion chamber 
as a result of mass transfer between cold and hot regions, in which 
any HC transferred are partially burned producing CO [23]. 

Apart from these inherent limitations, single-zone models can 
be used to provide an estimation of the ignition timing and the 
formation of NO, emissions. This is due to the fact that both of 
these phenomena are primarily controlled by the temperature of 
the hottest regions of the combustion chamber, which is close to 
the mean gas temperature of the mixture. Moreover, single-zone 
models can be used to assess the relative importance of chemical 
reactions and dominant species that are part of a chemical reaction 
mechanism describing the fuel-air oxidation. Thus, detailed reac- 
tion mechanisms can be compared to reduced ones under the same 
initial conditions and assumptions provided by the single-zone 
simulations. The latter features are especially important in the case 
of biofuels, the combustion chemistry of which is under continuous 
investigation. The relevant studies are perplexed by the variety of 
the biofuels, which include organic compounds of diverse homolo- 
gous series, such as alcohols, ethers, and esters. This is a significant 
departure from the combustion chemistry of the common pure 
hydrocarbons, which are usually used as surrogates for the simu- 
lation of diesel or gasoline combustion and do not include oxygen 


atoms in their molecules. Besides the study of the combustion 
chemistry of neat biofuels, single-zone models can be used to study 
the characteristics of various fuel blends or the effect of additives 
on biofuel HCCI combustion. This is an important feature, since the 
use of fuel blends and additives may provide the means for ignition 
timing and combustion rate control. Issues like the aforementioned 
provide the impetus for the ongoing HCCI research [88]. 


3.1. Hydrogen 


Goldsborough and Van Blarigan [89] studied numerically a free- 
piston hydrogen-fueled HCCI engine. The researchers modeled 
piston dynamics and the gas exchange process in their free-piston, 
two-stroke engine. They also used the HCT [90] chemical kinetics 
software to simulate a zero-dimensional thermodynamic process 
throughout the engine closed cycle. The hydrogen oxidation was 
accounted for using the chemical mechanism of Marinov et al. 
[67]. Heat transfer was accounted for by using Woschni's corre- 
lation [91], in which the heat transfer constants were modified 
to match experimental data from hydrogen combustion in a rapid 
compression machine. It was found that the engine operation was 
significantly dependent on the scavenging process, which altered 
the amount of fuel-air mixture trapped within the combustion 
chamber, the initial temperature of the charge, and the amount 
of fresh charge that was expelled directly through the exhaust port 
during this scavenging process. Overall, the engine was found to 
produce low NO, emissions and to provide almost constant volume 
combustion. The rapid combustion, which resulted in the simula- 
tions, could be partly attributed to the nature of the single-zone 
model used. 

Fiveland and Assanis [88] constructed a zero-dimensional 
numerical model for the simulation of the complete engine cycle 
using methane or hydrogen as a fuel. The hydrogen combustion 
chemistry scheme included 11 species and 25 chemical reactions, 
and was modeled using the CHEMKIN libraries [70] adapted for a 
variable volume plenum. Heat transfer was calculated using a sub- 
set of the k-e model under the isotropic turbulence assumption. 
This method yielded the characteristic velocity that was used to 
calculate the heat transfer coefficient. 

Simulation results for a hydrogen-fueled four-stroke engine 
with a compression ratio of 15 and 1.5 bar intake pressure required 
425K intake charge temperature to produce ignition near TDC. This 
high intake temperature requirement, at a moderate compression 
ratio and at a low fuel-air equivalence ratio of 0.3, was necessary 
since hydrogen ignition occurred at about 1050-1100 K in the sim- 
ulation. Under these conditions the researchers identified peaks in 
H, OH and HO; fractions near the ignition point. The simulation 
model was also used in a parametric investigation, in which it was 
found that decreasing the initial temperature or the fuel-air equiva- 
lence ratio, beyond a certain point, retarded ignition and eventually 
led to misfire. 


3.2. Alcohols 


Ng and Thomson [92] constructed a single-zone model to sim- 
ulate ethanol HCCI combustion and especially the effect of ethanol 
reformed gas on the ignition requirement and NO, emissions. They 
incorporated Marinov's set of chemical reactions [60] to simulate 
ethanol oxidation and integrated a NO, reaction mechanism from 
GRI-Mech [84] to simulate NO, formation. For the estimation of 
heat losses they used a convective heat transfer correlation (similar 
to Newton's Law of cooling). 

The researchers conducted a parametric study, in which the 
parameters varied were the fuel-air equivalence ratio, the EGR 
amount and the degree of reforming, i.e. the degree of ethanol con- 
version into CO and Hp. Thus, the initial charge contained ethanol, 
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CO, H2 and EGR at various proportions. It was found that fuel 
reforming resulted in an expansion of the HCCI operating region, 
the boundaries of which were determined by misfire and unac- 
ceptable NO, emissions. However, reforming tended to increase 
NO, emissions due to the resulting high peak combustion temper- 
ature. In relation to reformed gas, EGR was more effective toward 
expanding the simulated operating range at moderate equivalence 
ratios (0.2-0.5). 

Martinez-Frias et al. [93] investigated the direct use of wet 
ethanol using a simulation model. The use of wet ethanol as a fuel 
reduces the energy demands for its production, since the distilla- 
tion and/or dehydration processes, which are used to remove the 
water and provide neat ethanol, require a significant amount of 
energy (37% of the total energy in ethanol products for corn-based 
ethanol production). The objective of the study was to determine 
whether benefits could be obtained by the direct utilization of wet 
(as opposed to neat) ethanol into a turbocharged HCCI engine. The 
researchers constructed a simulation model for a TC HCCI engine 
fueled with wet ethanol (35 vol.% EtOH in H20). Apart from the 
TC, the full engine system model also incorporated an exhaust gas 
regenerator, a fuel vaporizer and a catalytic converter. 

The HCCI engine was modeled via a single-zone model, using 
CHEMKIN [94] to solve for chemical kinetics. The ethanol reac- 
tion mechanism was that of Marinov [60] and the NO,-related 
reactions were taken from GRI-Mech [84]. Heat transfer was also 
incorporated into the reciprocating engine model assuming a wall 
temperature of 430 K. Using this model, the researchers found that 
a minimum ethanol-in-water fraction of 35 vol.% (liq.) was neces- 
sary for a sustainable simulated HCCI operation. This provided an 
overall brake efficiency of 38.7% for the entire simulation system. 
The researchers estimated that this direct use of 35% wet ethanol 
significantly reduced the energy required for fuel production - rel- 
ative to neat ethanol production - by eliminating the need for 
dehydration and limiting the distillation needs to 32 of the total 
energy in ethanol products, in the case of corn-based ethanol pro- 
duction. 

In Ref. [95] Szybist used CHEMKIN-PRO [94] to simulate the HCCI 
combustion of n-heptane, isooctane, toluene and ethanol. The sim- 
ulation model was switched to the adiabatic mode. The mechanism 
used for the simulation of ethanol combustion was that of Mari- 
nov [60] including 56 species and 351 reactions. The author used 
this model to simulate HCCI combustion with the aforementioned 
fuels under normal naturally aspirated conditions with air dilution 
and under negative valve overlap (NVO) conditions. In the latter 
case, the mixture had a higher temperature at BDC (»500 K) and a 
significant amount of diluents such as O5, CO2, H20 and N2. These 
conditions resemble HCCI under NVO operation, in which a portion 
ofthe exhaust gases is trapped within the combustion chamber and 
raises the temperature of the incoming fuel-air charge. The results 
of the study showed that with low intake temperature (325 K) and 
no EGR, the most resistant fuel to ignite was ethanol requiring the 
highest compression ratio. However, under NVO conditions toluene 
was the most difficult to ignite. Moreover, with NVO the igni- 
tion timing between fuels converged to within a few CA degrees, 
with the exception of n-heptane at high O; dilution. On the con- 
trary, ethanol was essentially insensitive to the O5 dilution, under 
the conditions of the simulation, indicating a zero-order reaction 
rate with respect to O5. The researcher concluded that ignition 
in HCCI combustion under NVO conditions presented distinct fea- 
tures, which were not entirely consistent with the RON and MON 
parameters. 

Yang et al. [96] also investigated NVO HCCI combustion. They 
used a single-zone adiabatic model to simulate the reactions tak- 
ing place within the combustion chamber with the injection of 
gasoline, ethanol or methanol during the NVO period. The reac- 
tion kinetics was incorporated via the mechanisms of Marinov 


[60] for ethanol, and Lindstedt and Meyer [77] for methanol. The 
researchers studied both the thermal and the chemical effect of 
the fuel injection during the NVO period. They found that during 
this period both heat and active species such as Hz and CO can 
be produced, depending on the fuel type, oxygen availability and 
the in-cylinder temperature of the residuals during injection. At 
relatively high oxygen availability in the residuals the tendency 
was to produce more heat than active species, since the fuel was 
oxidized. In oxygen-poor residuals only part of the fuel was oxi- 
dized producing active species. Methanol and ethanol were less 
sensitive than gasoline to O; availability owing to their oxygenated 
nature. Methanol showed the greatest tendency in forming active 
species. 


3.3. Ethers 


In [48] Yao et al. used the SENKIN code of the CHEMKIN pack- 
age to simulate DME combustion. The single-zone model used 
was considered adiabatic and the chemical mechanism, which was 
used to describe the DME oxidation, was that of Lawrence Liver- 
more National Laboratory including 79 species and 336 reactions 
[52-55]. The scope of the study was to analyze the reaction path- 
ways of DME oxidation. The researchers validated the simulation 
results against experimental data, which indicated that the com- 
bustion rate and the peak combustion pressure were overestimated 
due to the nature of the single-zone model. However, the ignition 
timing was adequately predicted, implying a sufficient descrip- 
tion of the reaction pathways. It was found that DME presented a 
two stage auto-ignition with low temperature reactions (LTR) fol- 
lowed by high temperature reactions (HTR). HCHO and OH were 
the primary radicals during LTR. DME ignition was controlled by 
the H202 decomposition to two OH radicals, which provided a 
chain branching step. These hydroxyls were the primary attack- 
ers for the H atom abstraction from DME and the oxidation of 
CO. At relatively leaner DME concentration the LTR was not sig- 
nificantly affected but the CO oxidation deteriorated, indicating a 
temperature sensitivity of the CO oxidation process. This led to the 
generally accepted conclusion that CO emissions in HCCI engines 
originate from the relatively cold and/or lean regions of the com- 
bustion chamber, such as the crevices and the vicinity of the wall, 
in which the local temperature is adequately high to induce partial 
oxidation of the existing HC but not high enough to consume the 
resulting CO. 

Hofmann and Abraham [97] also used the internal combustion 
engine model included in the CHEMKIN package [94] to simulate 
n-heptane, DME and methyl decanoate HCCI combustion. The lat- 
ter was studied as a biodiesel surrogate. The simulation model 
incorporated Woschni's correlation [91] for the estimation of heat 
losses from the combustion chamber. The primary objective of the 
study was the investigation of ignition delay and NO formation 
for the three fuels. For this purpose, the researchers used Her- 
binet et al. [81] detailed chemical kinetic mechanism for methyl 
decanoate - a biodiesel surrogate - consisting of 3036 species 
and 8555 reactions, and Livermore National Laboratory scheme 
for DME oxidation including 79 species and 351 reactions [53-55]. 
NO formation was simulated using two alternative schemes; the 
extended reaction mechanism of Zeldovich [83] and the NO kinetics 
from GRI-Mech v.3 mechanism [84]. The Zeldovich NO formation 
scheme provided reduced NO concentration relative to the more 
extensive NO scheme of GRI-Mech for all fuels. Moreover, it was 
demonstrated that DME showed the shortest ignition delay fol- 
lowed by methyl decanoate and n-heptane. Thus, DME required 
higher engine speed and the lowest initial temperature. Although 
the primary consideration for NO formation was the peak combus- 
tion temperature, fuel chemistry was also found to play a secondary 
role. 
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3.4. Biodiesel 


Apart from Hofmann and Abraham mentioned in the preced- 
ing sub-section, biodiesel HCCI combustion was also simulated by 
Brakora and Reitz [98] using the SENKIN code of the CHEMKIN 
package [94]. The single-zone model was considered adiabatic for 
the purpose of the study, which was to investigate the formation 
of NO, from biodiesel- and diesel-fueled HCCI engines. N-heptane 
was used as the diesel surrogate, and a mixture of 2/3 n-heptane 
and 1/3 methyl butanoate was used to simulate biodiesel combus- 
tion. This proportion of n-heptane and methyl butanoate in the 
biodiesel surrogate was found to yield a better approximation for 
the actual C:H:O ratio found in the real biodiesel, relative to neat 
methyl butanoate or methyl decanoate. The latter two, being rela- 
tively short chain esters, provided a much higher C:O ratio than 
the actual biodiesel fuel. The chemical mechanism used for the 
biodiesel surrogate consisted of 56 species and 169 reactions [78]. 
NOx formation was incorporated via a 7 species and 19 reactions 
scheme, which was developed by Yoshikawa and Reitz [87] relying 
on reactions from the GRI-Mech mechanism [84]. 

Brakora and Reitz found that the ignition timing played a cru- 
cial role in the resulting NO, emissions. With an equivalent energy 
release and for the same ignition timing, the two surrogates pro- 
vided roughly the same peak combustion temperatures (within 
approximately 20K) and yielded similar NO, emissions at EVO. 
There was a difficulty in connecting the oxygen content of the fuel 
to the produced NO, emissions, especially at lean mixtures. Under 
stoichiometric conditions, the biodiesel surrogate yielded a 26% 
increase in NO, relative to the diesel surrogate. However, the pre- 
dicted difference in peak combustion temperature was about 13K 
lower for biodiesel, suggesting that the fuel-bound oxygen could 
play a role in the formation of NO, under limited O5 concentrations. 
It was also found that under the same energy content, ignition tim- 
ing and total N:O ratio in the mixture, for the two fuels, the diesel 
surrogate yielded higher NO, emissions. This suggested that hav- 
ing the same amount of O, available in the air, as opposed to being 
fuel-bound, increased the NO, emissions due to the alteration of 
the heat capacity of the mixture, which resulted in higher flame 
temperatures for the diesel surrogate. The study concluded that 
under HCCI conditions, in which near stoichiometric local condi- 
tions are not encountered, any NO, increase due to biodiesel would 
be small. 


3.5. Fuel blends and additives 


Apart from neat fuels, single-zone models have been used with 
fuels blends, to assess the effects of additives/blends on the HCCI 
combustion process. 

Szybist et al. [38] used a single-zone model to simulate the com- 
bustion of biodiesel blends in HCCI engines. In their study they used 
CHEMKIN [94] for the chemical kinetics simulation coupled with 
chemical mechanisms for n-heptane [99], which represented diesel 
fuel, and methyl butanoate [80], which accounted for biodiesel. The 
combined mechanism consisted of 670 species and over 3000 reac- 
tions, including NO, formation reactions, which were taken from 
GRI-Mech [84]. With this mechanism it was possible to study dif- 
ferent blends of n-heptane and methyl butanoate, which were used 
as surrogates for the simulation of diesel/biodiesel blends. Thus, the 
researchers simulated diesel/biodiesel blends of 0, 10 and 50% (vol.) 
in biodiesel with n-heptane/methyl butanoate blends, the relative 
proportion of which was chosen to match the oxygen content (wt%) 
of the experimental fuel blends. 

The simulation model predicted the effect of inlet temperature 
increase on the phasing and magnitude of the low temperature 
heat release (LTHR), and the phasing of the main combustion event. 
However, discrepancies were observed between experiments and 


simulation as regards the dependence ofthe results on the biodiesel 
concentration, and the duration and phasing of the LTHR and the 
main combustion event. These were attributed to the low cetane 
number of methyl butanoate, the approximation of the IVC tem- 
perature with the intake manifold temperature, and the inherent 
inability of the single-zone model to capture any fuel or tempera- 
ture stratification. The latter resulted to an overestimation of the 
peak combustion temperature and the combustion rate, which led 
to almost 100% combustion efficiency, very low CO and high NOx 
emissions. 

Saisirirat et al. [100] examined the effect of n-heptane blends 
with ethanol and 1-butanol, both experimentally and numerically. 
They used the SENKIN single-zone model of the CHEMKIN II code 
[94] to simulate n-heptane/ethanol and n-heptane/1-butanol fuel 
blends. The chemical mechanism used was that of Saisirirat et al. 
[44] consisting of 1046 species and 4398 reactions, which includes 
ethanol, 1-butanol and PRF oxidation sub-schemes. The simula- 
tion compared the HCCI combustion of neat n-heptane to the one 
obtained with the two fuel blends, i.e. n-heptane/ethanol and n- 
heptane/1-butanol both in 63:37 molar proportions. It was found 
that the addition of ethanol or 1-butanol to n-heptane affected the 
OH radical production and slowed down the chemical reactions 
during LTHR, decreasing the concentration of reactive intermediate 
species (HCHO, OH, HO;, etc.). Thus, alcohols induced a chemical 
effect on the LTHR. Owing to this reduction in chemically active 
species and to the low LTHR, which prevented a further tempera- 
ture rise of the mixture, the main combustion event was delayed 
with the addition of alcohols. 

Mack et al. [101] constructed a single-zone model for the sim- 
ulation of ethanol/diethyl ether (DEE) blends. The single-zone 
model incorporated Woschni's correlation for the estimation of 
heat losses; it used a DME-derived chemical mechanism [52] for 
the chemistry of DEE oxidation, and Marinov's ethanol oxidation 
scheme [60] for the combustion chemistry of ethanol. The com- 
bined mechanism consisted of 112 species and 484 reactions. The 
purpose ofthe study was to investigate the relative reactivity of DEE 
and ethanol in a fuel blend of 252; DEE in ethanol. The experiments 
conducted with carbon-14 isotope tracing showed that DEE was 
more reactive than ethanol, yielding more CO; relative to ethanol. 
The simulation results supported this conclusion, since it was found 
that DEE began to react earlier than ethanol. However, a signifi- 
cant elongation of the heat release, which would be desirable to 
control the combustion process, was not observed using the fuel 
blend. 

The same research group also studied the effect of small 
amounts of di-tertiary butyl peroxide (DTBP) (up to 3%) on the HCCI 
combustion of ethanol and diethyl ether blends [45]. The intention 
wasto elongate the combustion event using the exothermic decom- 
position of DTBP to promote early DEE ignition. The single-zone 
model used was the one mentioned earlier and the combustion 
reaction scheme was enhanced by the adiabatic decomposition 
of DTBP based on the work of lizuka et al. [102]. It was found 
via simulation, and verified by experiments, that the addition of 
DTBP advanced the ignition timing in the DEE-ethanol blend. The 
simulations showed that the decomposition of DTBP raised the 
temperature ofthe mixture and advanced ignition, although poten- 
tial chemical effects of DTBP were not excluded. 

Yao et al. [103] used a zero-dimensional thermodynamic model 
coupled with CHEMKIN III to simulate a DME/methanol blend in an 
HCCI engine. The chemical mechanisms used were the ones of Cur- 
ran et al. for DME [52], consisting of 78 species and 336 reactions, 
and Held and Dryer [74] for the reactions of methanol, consisting of 
22 species and 89 reactions. The authors compared the simulation 
predictions to experimental results of pressure and operating range 
in Apme — Ameou map. The analysis ofthe simulation results showed 
thatin DME/methanol dual fuel oxidation the low temperature heat 
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release of DME was inhibited, resulting only in a high temperature 
heat release. The OH radicals were mainly formed by the decompo- 
sition of H203 at blue-flame reactions, and by H * O? 2» O* OH and 
O*H50- OH * OH at hot flame reactions. 

Similar results were found by the same Institute in another study 
[104]. The simulation model used has already been described in the 
first part of their paper [48] presented in Section 3.3. In Ref. [104] 
the researchers studied the HCCI combustion of a DME/methane 
blend. Under the simulation conditions, the dual fuel blend showed 
a two stage autoignition scheme. This was different from the clas- 
sical two stage ignition observed in neat DME oxidation, since it 
corresponded to two high temperature heat release rates, one of 
DME, which was found to oxidate first, and the other of methane. 
It was observed that methane - like methanol - inhibited the LTHR 
but promoted the high temperature heat release (HTHR) of DME. 
On the other hand, DME promoted methane oxidation both due to 
thermal and chemical considerations. The thermal effect was due 
to the temperature increase, which was the result of DME oxida- 
tion, while the chemical effect was due to HCHO and H202, which 
were formed during DME oxidation and increased the OH radical 
pool. The researchers also identified the major reaction paths for 
the dual-fuel oxidation. 

In [105] Ogawa et al. studied the suppression effects of addi- 
tives such as water, methanol, ethanol, hydrogen and methane, on 
the combustion of a DME fueled HCCI engine. The researchers used 
SENKIN under adiabatic conditions for the single-zone simulations 
and the corresponding chemical mechanisms were adapted from 
Curran et al. [52] for DME, Marinov [60] for Ethanol, and Westbrook 
etal. [75,76] for Methanol. In the experiments the suppressors were 
injected directly into the combustion chamber, while the DME was 
introduced into the air stream at the intake pipe. The direct injec- 
tion of suppressors induced a cooling effect on the charge. Water 
injection had a predominantly cooling effect, i.e. it lowered the tem- 
perature of the mixture during its vaporization, thereby retarding 
ignition. Simulations conducted with the assumption that gaseous 
water was introduced in the mixture at BDC, resulted in a slight 
delay in radical formation and suggested that the water addition 
also had a minor chemical effect. Similar findings characterized 
methanol addition. Its main effects were a thermal one observed 
during the experiments and a chemical one suggested by simula- 
tions. 

As in the case of water addition simulations, methanol addi- 
tion simulations were conducted assuming that methanol existed 
in gaseous state at the beginning of calculations. It was found that 
methanol addition severely suppressed LTHR to near-zero values, 
significantly reduced OH, H and O radicals, and retarded main com- 
bustion. Further simulations with ethanol, hydrogen and methane 
addition showed that LTHR was slightly reduced with hydrogen or 
methane addition, and essentially disappeared with the addition 
of ethanol or methanol, in order of significance. Moreover, it was 
found that despite the different simulation conditions with respect 
to additive and its quantity, the OH concentration at the onset of 
LTHR remained fairly constant. This suggested that the chemical 
effect induced by the additives was linked to a reduction in the OH 
radicals prior to LTHR. 

This assumption regarding OH radicals was supported by later 
studies with Hj and MeOH addition to DME-air mixture. Specifi- 
cally, Shudo et al. [106,107] investigated the effect of H5 addition, 
which originated from on-board methanol reformed gas, on the 
DME HCCI combustion. The concept of the combustion engine sys- 
tem and the experimental setup are given in Fig. 6(a) and (b). 

The SENKIN [94] simulation model was used to describe the 
effect of hydrogen addition, and the reaction mechanism used was 
that of Fischer et al. [53-55] consisting of 79 species and 351 reac- 
tions. The simulation results showed that the addition of H3 in the 
DME-air mixture suppressed DME consumption. This was due to an 
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Fig. 6. (a) Concept of HCCI combustion engine system fueled with DME and MRG 
onboard-reformed from methanol utilizing engine exhaust gas heat [107] and (b) 
experimental setup [106]. 


antagonistic effect for the consumption of OH radicals; in the pres- 
ence of H5, OH radicals reacted preferentially with H5 instead of 
DME, to yield the relatively stable — at low temperatures - interme- 
diates HO; and H50;. This led to a suppression of LTHR, a decrease 
in the in-cylinder temperature rise and a consequent retardation of 
the main combustion event. This thermal effect was the main rea- 
son for the combustion retardation, since the ignition temperature 
was not altered appreciably. These findings implied that hydrogen 
or hydrogen-containing additives, such as methanol reformed gas 
(MRG), could be used to control the ignition timing of DME. The 
researchers also studied the conversion of methanol to DME and H3, 
via dehydration- and thermal decomposition-methanol reformers, 
respectively, to conclude that the H2/DME fraction in MRG could be 
continuously controlled by adjusting the methanol supply in each 
reformer. The addition of CO in neat DME, which was studied in 
[106], was also found to suppress the DME oxidation at low tem- 
perature, but its effect was less pronounced than the H3 addition 
effect. 

Yamada et al. [108] studied the effect of methanol or O3 addi- 
tion to DME-air mixture under HCCI combustion conditions. They 
used the SENKIN [94] code to simulate the engine operation, and 
implemented a chemical mechanism by Curran et al. [52] to account 
for the DME and methanol addition chemistry. Ozone addition was 
accounted for by including the thermal decomposition reaction of 
ozone according to Heirmel and Coffee [109]. The experimental 
and simulation results showed that methanol addition retarded the 
main heat release event significantly; 8% addition of methanol in 
DME resulted in a 15° CA delay. It was suggested that this delay was 
due to the competitive consumption of OH by methanol instead of 
DME, which led to the formation of relatively stable intermediates 
such as HCHO and HO;. In contrast, the addition of low O3 amounts 
(0.015% in DME) advanced the ignition timing by as much as 20° CA. 
This was due to the decomposition of O3 at relatively low temper- 
atures (—600 K), which provided oxygen atoms and activated the 
mixture. This supply of oxygen atoms at temperatures lower than 
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Fig. 7. The First Law of Thermodynamics applied to a zone of a multi-zone model. 


those of normal cool flames suppressed the production of formalde- 
hyde, which was a chain termination step under low temperatures. 
Thus, more heat was liberated during cool flame reactions, leading 
to an earlier thermal ignition. 


4. Multi-zone modeling 


The limitations encountered with single-zone models can be 
partly overcome with the use of multi-zone models. The latter 
assume that the combustion chamber is divided into a sufficient 
number of regions or zones, each possessing uniform thermody- 
namic properties. The sufficiency of the number of zones is usually 
determined by a sensitivity analysis, i.e. the minimum number of 
zones used is defined as the number beyond which the simulation 
results tend to converge. 

The temperature, composition, volume and all thermodynamic 
properties of each zone are considered uniform within the bound- 
aries of the zone, but can be different from the thermodynamic 
properties of other zones. This provides an opportunity to study 
the effects of temperature or species stratification within the com- 
bustion chamber, and usually leads to a better estimation of the 
combustion duration, the pressure rise rate and the emissions for- 
mation. Mass and heat transfer within the combustion chamber 
play a critical role in the estimation of the aforementioned operat- 
ing characteristics and since the combustion chamber is considered 
an ensemble of thermodynamically distinct zones, the processes 
of mass and heat transfer between zones must also be described 
by the multi-zone model. Fig. 7 shows the first law of thermo- 
dynamics applied to a zone, which in the general case includes 
heat and mass transfer. As regards mass transfer, each zone can be 
thermodynamically defined as a closed or an open thermodynamic 
system, depending on whether mass flow between zones is pro- 
hibited or allowed, respectively. Likewise, heat transfer between 
zones and to the combustion chamber wall can also be permitted 
or neglected. 

However, the assumptions regarding mass and heat transfer 
affect the results of the simulation significantly; heat transfer 
affects the temperature stratification within the combustion cham- 
ber, the colder regions of which are sources of unburned HC and CO. 
The formation of these pollutants is also affected by mass trans- 
fer between zones, since any HC or CO flowing from a colder to a 
hotter region/zone is expected to be (partially) oxidized [110]. The 
predicted combustion rate, pressure rise rate, combustion duration 
and peak pressure are also affected by the choice of the sub- 
models used to describe heat and mass transfer. The multi-zone 
models examined subsequently are divided into stochastic based 
multi-zone models, in which some properties of the mixture are 
considered random variables, and phenomenological multi-zone 
models, in which probability considerations have been excluded. 


In [62] Mosbach et al. constructed a stochastic reactor model 
(SRM) based on a probability density function (PDF) approach. The 
PDF was approximated by an ensemble of 100 stochastic particles. 
These did not contain any spatial or geometric information. The 
random variables were the species concentration and the tempera- 
ture. The model included heat transfer to the combustion chamber 
wall using Woschni's correlation [91] and the magnitude of the 
temperature fluctuations was regulated by a coefficient. Turbulent 
mixing was accounted for via a mixing intensity term that included 
the characteristic mixing time. The fuel simulated was a blend of 
EtOH/DEE, although a validation of the model was also presented 
with neat EtOH as fuel. The validation showed good agreement 
between experimental and simulated results as regards pressure 
traces. The model also captured qualitatively the experimental vari- 
ation of CO and HC emissions with the combustion phasing, i.e. 
the angle at which 50% of the total heat has been released. In 
these simulations, the equivalence ratio remained constant and 
the combustion phasing was regulated via the initial temperature. 
The stochastic reactor model was also used to study numerically 
the effect of the DEE in EtOH ratio on the combustion phasing, 
the combustion duration and the emissions of the HCCI engine. It 
was shown that as the fraction of DEE in the DEE/EtOH blend was 
increased, the combustion phasing advanced before TDC, the com- 
bustion duration was shorter, and the CO emissions were reduced. 
This fact provided indication that DEE could be used with ethanol 
to regulate the combustion phasing and duration of the HCCI com- 
bustion. 

Noda and Foster [111] numerically studied hydrogen HCCI com- 
bustion by constructing a multi-zone model; the different zones 
catered for any temperature or composition inhomogeneities. The 
zones were not assigned a specific location within the combus- 
tion chamber, and mass exchange between them was neglected. 
The heat losses were calculated via the Woschni [91] correlation 
using the combustion chamber area and the mean gas tempera- 
ture. The combustion chemistry was accounted for by incorporating 
11 species and 23 reactions from Westbrook and Dryer [69] for 
the oxidation of hydrogen, and the extended Zeldovich mecha- 
nism of Lavoie et al. [83] for NO, formation. The researchers used 
10-30 zones to examine the effect of temperature and fuel inho- 
mogeneities on HCCI combustion. Both temperature and - to a 
much lesser degree - fuel inhomogeneities increased the com- 
bustion duration compared to single-zone simulations. The mean 
gas temperature and the type of temperature inhomogeneity also 
affected combustion duration and the formation of NO emissions. 
Provided the mean gas temperature was high enough to avoid 
bulk quenching of the mixture, the most favorable type of tem- 
perature stratification included a small amount of hot gases to 
induce autoignition and avoid excessive NO formation and a greater 
amount of colder mixture. This type of stratification expanded the 
operating range in a mean gas temperature versus temperature 
inhomogeneity map. 

Komninos et al. [112] and Komninos and Rakopoulos [40,113] 
used a multi-zone model to simulate hydrogen [112] and ethanol 
[40,113] HCCI combustion. In all studies, the multi-zone model 
used divided the combustion chamber into several zones, which 
occupied specific geometric regions (Fig. 8). 

The only uniformity assumed was that of pressure, while heat 
and mass transfer between zones was accounted for, taking into 
account the geometric configuration of the zones; each zone could 
communicate with its immediate inner and outer zone (Fig. 9). The 
model also accounted for a crevice zone, which was in direct com- 
munication with the outmost zone of the combustion chamber. In 
Ref. [112], 11 zones were used - including the crevices - to simulate 
hydrogen HCCI combustion and the Annand correlation [114] was 
employed to estimate the heat flux from the outmost zone to the 
combustion chamber wall. The chemistry of hydrogen oxidation 
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Fig. 8. Zone configuration in the multi-zone model [112]. 


was accounted for via the mechanism of Marinov et al. [54] involv- 
ing 10 species and 19 reactions, augmented by a NO, formation 
scheme adapted from the isooctane oxidation reaction mechanism 
of Golovitchev [86]. A parametric investigation was conducted with 
the multi-zone model to determine the effect of various operating 
parameters on hydrogen HCCI combustion. Compression ratio and 
initial temperature were the dominant factors affecting the igni- 
tion timing and pressure rise rate. Engine speed also affected the 
ignition timing, mainly due to the variation of the available time for 
the occurrence of chemical reactions. It was found that increasing 
the load, i.e. the equivalence ratio, NO, emissions were increased 
significantly. Inlet pressure increased almost proportional to the 
attainable imep, without compromising the fuel conversion effi- 
ciency for the conditions simulated. 

In [40] Komninos and Rakopoulos numerically compared 
ethanol and isooctane HCCI combustion. The multi-zone model 
used was the one described in the preceding paragraph with some 
modifications; the zones number was increased to 16, the zone 
resolution was refined near the combustion chamber wall, and the 
wall heat flux estimation approach abandoned the Annand corre- 
lation and was based on the temperature gradient near the wall. 
These modifications were necessary to better estimate the forma- 
tion of HC and CO emissions [115]. The chemical mechanism used 
for the ethanol oxidation was that of Róhl and Peters [57] consist- 
ing of 38 species and 228 reactions. Since this mechanism did not 
include NO, reactions, it was augmented by 13 reactions involving 
nitrogen oxides formations [86]. The model vas validated against 


Fig. 9. Schematic of heat transfer between zones and to the cylinder wall [113]. 


experimental results and captured adequately the pressure traces 
and HRR for both ethanol and isooctane. 

According to the simulation results, CO emissions were formed 
near the combustion chamber walls, thereby justifying the finer 
zonal resolution used in these regions (see Fig. 10). As shown in 
Fig. 10, the amount and location of CO production was found to 
depend on load; at higher loads, CO production advanced to the 
inner zones. 

Total CO emissions were adequately captured but HC emissions 
were overestimated, especially for ethanol. An investigation for the 
composition of the HC emissions for both isooctane- and ethanol- 
fueled HCCI combustion showed that in the ethanol cases examined 
the vast majority of unburned HC consisted of oxygenated HC 
species (OHC). These were primarily ethanol, acetaldehyde and 
formaldehyde. The predicted percentage of OHC in total HC (aprox. 
95% by mass) for ethanol was not very sensitive to the load 
variation under the simulation conditions. On the contrary, the 
relative amount of OHC in total HC was lower (19-33% by mass) 
in the case of isooctane and decreased with the increase of load. 
It was concluded and demonstrated that this would induce a 
significant error in the HC measurement if a measuring device 
such as the FID was used, since the FID principle possesses an 
inherent inability to accurately detect oxygenated carbon atoms 
[116]. 

Using the same model in [113] the authors compared two 
ethanol oxidation mechanisms, namely Marinov's [60] and Róhl 
and Peters' [57]. Marinov's mechanism consisted of 57 species and 
383 reactions, while Róhl and Peters' mechanism, which was a 
reduced version of Marinov's, consisted of 38 species and 228 reac- 
tions. Both mechanisms were augmented by NO, related reactions 
taken from [86]. The results showed that the two mechanisms 
predicted essentially the same dominant HC species, with the 
exception of formic acid (HCOOH) and to a lesser extent methanol 
(CH30H), both of which were not included in the reduced mech- 
anism. The mechanism reduction had a significant impact on 
the computational time: although the species used in Róhl and 
Peters' mechanism were 2/3 of the species used in the Marinov 
mechanism, the Röhl and Peters’ mechanism required half the com- 
putational time. The simulation results showed a slightly higher 
predicted combustion rate by the Róhl and Peters' mechanism rel- 
ative to Marinov's mechanism, despite the fact that the predicted 
initial heat release was the same for the two mechanisms. The 
resulting temperature fields were almost identical as shown in 
Fig. 11. Therefore, the two chemical mechanisms did not induce 
significant differences in the predicted performance. NOx, CO and 
HC emissions were essentially the same for the two mechanisms, 
with the greatest disagreement occurring between predicted CO 
emissions (~9%) at the low load cases. This was connected with the 
slightly lower combustion rate predicted by Marinov. 

Sjoberg and Dec [117,118] used a multi-zone model to study the 
effects of engine speed, equivalence ratio, initial temperature and 
boost pressure on ethanol HCCI combustion. They used the multi- 
zone model provided by the Senkin application of the CHEMKIN 
package [94], which treated each zone as a single lumped mass 
with uniform composition. The only interaction between zones 
was through pressure, which was considered uniform. Thus, the 
volume of each zone depended on its mass and its temperature 
relative to the other zones. The 11 zones used were adiabatic in 
principle. This led to the need to introduce initial temperature 
stratification for the 10 active zones to match the experimental 
pressure trace. The 11th zone was assigned a temperature low 
enough to render it chemically inactive throughout the engine 
cycle. This zone included 16% of the total mass and compensated 
for the combustion inefficiencies, heat transfer and blow-by, which 
were not directly modeled. The simulated compression ratio was 
also reduced relative to the experimental one to account for the 
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Fig. 10. CO simulation results for ethanol; left: in-cylinder CO emissions and right: zone CO production rate [40]. 


neglected heat losses. The thermal stratification of the zones was 
supported by KIVA simulations and reproduced adequately the 
variation of combustion phasing with temperature. The chemical 
mechanism used consisted of 58 species and 310 reactions [61], 
and was part of a larger gasoline surrogate mechanism produced 
by the joint effort of the National University of Ireland-Galway 
and Lawrence Livermore National Laboratories [119]. 

The simulations above showed that the ethanol oxidation mech- 
anism reproduced adequately the effect of equivalence ratio, engine 
speed and boost pressure on ignition. The same model was used by 
these researchers to study the effect of EGR and its constituents on 
the HCCI autoignition of ethanol [118]. In this study, it was found 
that the addition of simulated EGR did not affect significantly the 
ethanol autoignition timing for the conditions examined. The effect 
of EGR on other fuels, such as isooctane and PRF, was more pro- 
nounced. This was mainly owed to the low sensitivity of ethanol 
oxidation on oxygen concentration, which was reduced with the 


increase of simulated EGR. Moreover, it was suggested that the 
presence of water in the EGR promoted ethanol auto-ignition. 


5. Multi-dimensional modeling 


Besides the aforementioned favorable attributes of the multi- 
zone models, they also have limitations. Even in the cases where 
spatial considerations have been accounted for by multi-zone mod- 
els, they usually do not allow high local resolution. Moreover, 
there is an inherent difficulty to fundamentally describe the motion 
of the charge within the combustion chamber and the effects of 
turbulence. The problem is even more pronounced in the case 
of Diesel-like HCCI engines, in which the mixture formation is 
achieved via direct fuel injection within the combustion cham- 
ber. Direct fuel injection imposes new demands on any simulation 
model, since the mixture - and the resulting temperature - strati- 
fication within the combustion chamber is affected significantly by 
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Fig. 11. Zone temperature distributions produced by Róhl and Peters' (solid line) 
and Marinov's (dashed line) mechanisms for the high load ethanol case ($ = 0.2641) 
[113]. 


the injection characteristics. Phenomena such as wall wetting, fuel 
dispersion, fuel vaporization under the low pressure-low temper- 
ature conditions encountered during early direct injection, must 
be described. In these cases, the incorporation of accurate models 
for the estimation of the physicochemical fuel properties becomes 
crucial, especially for biofuels, since their development and use has 
a relatively short history. 

The aforementioned issues can be dealt with by using multi- 
dimensional models, i.e. CFD models, which describe in detail the 
charge motion within the combustion chamber and the turbulence 
effects, provide a superior spatial resolution of the combustion 
chamber, and can account for the description of DI HCCI combustion 
by incorporating the dynamics and kinematics of the fuel injection. 

In the remaining part of this section, direct injection HCCI CFD 
models will be discussed first while dedicated HCCI CFD models 
assuming homogeneous conditions at IVC will be presented subse- 
quently. 


5.1. Direct injection HCCI 


Kim et al. [120] investigated the influence of the spray angle 
and advanced injection timing on the combustion performance 
and emission characteristics in a DME fueled compression ignition 
engine. They used the KIVA software [121], which was supple- 
mented with liquefied DME fuel properties and their temperature 
dependence [122-126]. The oxidation process of the DME fuel 
was described with a detailed chemical kinetic model [53-55], 
which consisted of 351 reversible elementary reactions among 79 
species. This mechanism was expanded by including NO, forma- 
tion reactions taken from a modified version of GRI Mech [84]. A 
two-step phenomenological soot model [127] was also applied to 
calculate soot emissions. The DME spray and atomization char- 
acteristics were simulated via the Kelvin-Helmholtz (KH) and 
Rayleigh-Taylor (RT) hybrid breakup model [128]. Turbulence was 
accounted for via the modified Renormalization Group (RNG) k-& 
model [129], and wall impingement was simulated via the spray 
wall interaction model [130] including droplet splash/spreading. A 
computational grid with a unit cell size of about 2 mm? was formed 
for the evaluation of spatial distribution in the cylinder and the 
computational grid was divided into a 60° sector mesh, due to the 
symmetric arrangement of the six-hole injector. The calculation 
time step was set at 0.1 ws. 
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Fig. 12. Cylinder locations according to the injection timing and five injection angles 
[120]. 


The model was validated against experimental measurements 
at three injection timings, with an injection angle of 156*. The 
experimental and simulation results at an injection timing of 60? 
BTDC had two peaks in the heat release rate, indicating low- and 
high-temperature DME oxidation. The validated model was subse- 
quently used for a numerical study aimed at analyzing the effect 
of the injection angles and timings shown in Fig. 12. The result- 
ing fuel spray patterns created by these configurations are shown 
in Fig. 13. These patterns revealed fuel spray-wall interactions and 
accumulation of fuel in the squish area and, possibly, the crevice 
regions. 

Generally, the fuel injected vaporized rapidly as it entered the 
combustion chamber due to the high volatility of DME. At an early 
injection timing (60° bTDC), the variation of the injection angle 
affected significantly the pressure and heat release rate diagrams; 
at wider injection angles the peak pressure and HRR were both 
reduced, since some of the injected fuel reached the relatively cold 
crevice and combustion chamber wall regions. At later injection 
timings, the effect ofthe injection angle was much less pronounced, 
since the fuel was injected within the piston bowl. As the injection 
timing approached TDC, the combustion mode shifted progres- 
sively from partially homogeneous to diffusive combustion. Soot 
emissions were very low for all cases examined due to the chemi- 
cal structure ofthe fuel, which excludes C-C bonds. Another reason 
was the relatively premixed nature of the mixture, especially at 
early injection timings. 

NO, emissions were very low for the early injection timings (IT 
40° bTDC and earlier) due to the premixed nature of the mixture. 
As the injection timing approached TDC, NO, emissions increased, 
probably owing to the predominance of the in-cylinder inhomo- 
geneities. At advanced injection timings and/or wide injection 
angles, the HC and CO emissions were higher due to the introduc- 
tion of fuel into the squish and crevice regions. At narrow angles, the 
high HC and CO emissions were due to the fuel impingement on the 
piston bowl. As regards performance, a narrow angle produced low 
ISFC at relatively advanced injection timing. Other studies focusing 
solely on the spray characteristics of DME by the research group of 
the Hanyang University can be found in [131,132] 

Chen et al. [133] investigated experimentally and numerically 
the effect of methanol injection timing in a DME/MeOH dual 
fuel engine, in which the DME was introduced at the inlet port. 
For the numerical simulation of the closed part of the engine 
operation, the researchers coupled chemical kinetics with the CFD 
code FIRE [134]. The chemical mechanism used to estimate the 
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Fig. 13. Effect of the injection angle and timing on calculated spray and evaporation 
distributions in the cylinder according to the crank angle [120]. 


methanol/DME reaction rates consisted of 30 species and 44 reac- 
tions including the fuel chemistry [50], the Zeldovich mechanism, 
and the N50 intermediate mechanism [82]. The effect of turbulence 
was accounted for by adopting a compressible version of the stan- 
dard two-equation k-e turbulence model. The source terms in the 
transport equations of FIRE were calculated with the assumption 
thatthe reaction rate is mainly determined by kinetic and turbulent 
time scales. Major submodels used in the computation included 
the WALLJET1 wall interaction model, the Dukowicz evaporation 
model, and the HuhGosman breakup model [135]. The Han-Reitz 
wall heat transfer model was used to estimate the wall heat flux 
[136]. The spray angle was 152° formed by a 7-hole injector, which 
allowed for the simulation of a 51? sector consisting of 30,000 
grid cells at BDC. The simulations examined the effect of MeOH 
injection at two injection timings, namely —26? and —6° aTDC. 

It was found that methanol injected at —26° aTDC had no 
effect on the low temperature reaction of DME but accelerated the 
heat release of the high temperature reaction, producing locally 
high temperature regions and increasing NO concentration. In this 
case, the injected methanol was completely oxidized. At the later 
injection timing (—6° aTDC) the DME found in the crevices was 
higher, and the high temperature region was dispersed in the 
clearance volume and the region near the combustion chamber 
wall, leading to lower NO concentration. Some methanol vapor 
remained in the combustion chamber prolonging the combustion 
duration. This led to an investigation aimed at decreasing the com- 
bustion duration at this - relatively late — injection timing. The 
researchers studied numerically the effect of injection pressure on 
the combustion duration, while maintaining the injection timing 


a 


Fig. 14. Computational mesh of combustion chamber at TDC used in [138]. 


at —6° aTDC. They found that increasing the injection pressure from 
20 to 30MPa shortened the combustion duration, increased the 
thermal efficiency, and induced an increase in peak temperature 
and NO concentration. However, it was difficult to experimentally 
examine the effect of increased injection pressure due to the low 
viscosity of methanol. 


5.2. Fully premixed HCCI 


In Ref. [137] Yao et al. examined the HCCI combustion of a dual 
fuel engine fueled with DME/MeOH blends. They used the STAR- 
CD/KINETICS CFD code to simulate the combustion of the charge 
that was uniform at IVC. The RNG k-e model was included to 
account for the effects of turbulence, and the PISO algorithm was 
used for the transient flow of the engine. The CFD code was coupled 
with the CHEMKIN code to obtain the species concentration evolu- 
tion, and the DME/methanol oxidation was described by a reduced 
mechanism consisting of 27 species and 35 reactions [50]. The time 
step corresponded to 0.1? CA at an engine speed of 1400 rpm. The 
researchers validated the model results against an experimentally 
obtained pressure diagram and found adequate agreement with 
the measurements. Under the same experimental conditions they 
also conducted an investigation into the formation of HC and CO 
emissions. They found that UHC emissions mainly consisted of the 
unburned fuels and CH20. The unburned DME and methanol mainly 
resided in the piston-ring crevice region, while most of the CH20 
emissions were located next to the cylinder-liner wall. CO emis- 
sions were mainly found in the regions near the top surface of the 
piston. 

The authors also studied the effect of load (in terms of equiva- 
lence ratio) and the relative proportion of DME in the fuel blend on 
the HC and CO emissions. Increasing the DME ratio, while main- 
taining the overall equivalence ratio, improved the combustion 
efficiency and reduced the HC and CO emissions. Moreover at all 
loads examined, a critical value for the DME ratio was found; when 
reducing the ratio below this value induced a significant increase in 
both CO and HC due to incomplete combustion. However, the criti- 
cal DME ratio was not the same for all equivalence ratios examined. 
At high equivalence ratios, the critical DME ratio shifted toward 
higher values due to the increase in methanol concentration at 
high loads, which induced an inhibiting effect on DME ignition. 
A critical value was also found for the peak combustion temper- 
ature (1400 K), above which the HC and CO emissions were greatly 
reduced. When the load, i.e. the overall equivalence ratio, was 
increased while maintaining the DME proportion, HC and CO emis- 
sions were reduced. 

Using roughly the same simulation model, Huang et al. [138] 
studied numerically the HCCI combustion of a DME fueled engine. 
A reduced chemical mechanism for DME HCCI combustion pro- 
cess was adopted, which was proposed by Yao et al. [48] and Liang 
et al. [49]. The reduced chemical mechanism was derived from the 
detailed chemical mechanism and consisted of 26 species and 28 
reactions. A single layer grid was applied with an average cell size 
of about 2 mm, yielding about 1658 cells at IVC in the combustion 
chamber (Fig. 14). 


N.P. Komninos, C.D. Rakopoulos / Renewable and Sustainable Energy Reviews 16 (2012) 1588-1610 1605 


Initial locations 
for LTR 


Initial locations 
for HTR 


650K 1650K 


Fig. 15. Initial locations for LTR and HTR in [138]. 


The computational time for the closed part of the engine cycle 
was approximately 18h in a Pentium IV 3.0 GHz CPU PC with 1 GB 
memory. Under homogeneous charge IVC conditions, the DME HCCI 
combustion process does not take place in the entire combustion 
chamber simultaneously. The low temperature reaction was initi- 
ated primarily at the region near the piston surface and the squish 
region, and subsequently moved toward the cylinder head. The 
high temperature reaction was initiated in the combustion cham- 
ber core zone and squish region, and next proceeded to the rest of 
the combustion chamber (Fig. 15). 

The majority of HC emission was unburned fuel and CH20, 
and resided with CO in the bottom, middle, and upper part of 
the piston-ring crevice regions, respectively. At very low loads, 
DME was partially oxidized producing unburned HC and CO. As 
the load increased, the unburned HC decreased and the CO emis- 
sions increased and then decreased. This behavior can be explained 
considering the successive steps through which the hydrocarbon 
oxidation takes place, i.e. the oxidation follows the general scheme 
HC CO — CO;. At extremely low loads the fuel and any unburned 
HC produced by its decomposition fails to be converted to CO to a 
great extent. At higher loads, the HC to CO conversion takes place 
more vividly and, at the highest attainable loads, both HC and CO are 
oxidized to produce CO; through the aforementioned generalized 
oxidation scheme. 

Chun-lan et al. [56] simulated neat DME HCCI combustion. They 
used a single-zone model to compare two chemical mechanisms: 
a detailed mechanism of 97 species and 457 reactions, which was 
obtained by combining the mechanism from Curran et al. [54] with 
the NO, formation reactions of Marinov [85], and a simplified reac- 
tion mechanism of 36 species and 73 reactions. The latter included 
alow temperature, a negative temperature coefficient, and a pyrol- 
ysis and oxidation submodel along with the reaction path for NOx 
formation. 

The FLUENT CFD software was also used to obtain insight into 
the formaldehyde and formic acid formation process during the 
DME HCCI combustion. For this reason the CFD code was coupled 
with the three submodels, i.e. the low temperature, the negative 
temperature coefficient, and the pyrolysis and oxidation submod- 
els. The RNG k-e turbulence model and the EDC combustion 
model were implemented in the CFD code. Using this model, the 
researchers found that both formaldehyde and formic acid existed 
in the exhaust. Their concentration was close to the experimen- 
tally determined, in the case of formaldehyde, but the measured 
formic acid was much less than the one predicted. The CFD cal- 
culations showed that the formation of these species was located 
at the region near the combustion chamber wall and above the 
ring crevices. Their formation originated from the partial oxida- 
tion of DME, which had entered the crevices during combustion 
and flowed back to the combustion chamber during expansion. The 
partial oxidation of DME gave rise to formaldehyde and formic acid, 
the further oxidation of which was suppressed due to the low local 
temperatures induced by the gas expansion and heat transfer to the 
cylinder wall. Some of the formaldehyde and formic acid formed 
were transferred to the central part of the combustion chamber, 


due to mass transfer, indicating that their concentration can be 
affected by the combustion chamber configuration and the exist- 
ing flow field. At EVO a relatively high concentration of these gases 
remained near the wall. 

In Ref. [139] Park numerically analyzed the combustion and 
emission characteristics and the optimal operating conditions of 
homogeneous charge compression ignition (HCCI) engines fueled 
with dimethyl ether (DME). For this investigation, a modified 
version of KIVA-3V code [121] was used to simulate DME HCCI 
combustion with a two-dimensional computational mesh. The 
CHEMKIN chemistry solver [70] was integrated into KIVA-3V, in 
order to solve the chemistry and the chemical mechanism used for 
the simulation of DME consisted of 351 reactions and 79 species 
[53-55]. A two-step phenomenological model was used to simu- 
late soot emission, since the study extended the investigation into 
rich DME combustion. The soot model considered both soot for- 
mation, using a modified Hiroyasu model that considers acetylene 
(C2H2) as a soot precursor [140], and soot oxidation employing 
the Nagle-Strickland-Constable model [141]. For NO, emissions, 
a mechanism derived from the Gas Research Institute (GRI) NO 
mechanism [84] was used. 

During the simulations, the charge was assumed homogeneous 
at intake valve closure. The researcher generated an equivalence 
ratio versus peak temperature operating map based on numeri- 
cal simulations, and compared the DME operation to n-heptane 
operation. High CO concentration was found at low peak combus- 
tion temperatures and rich mixtures. CO emissions for DME were 
greater than those for n-heptane, because of a lower initial mole 
fraction of oxygen. HC emissions for DME and n-heptane were quite 
similar at lean mixtures. However, HC emissions for DME were 
lower than those for n-heptane below a peak cycle temperature 
of 1600 K and an equivalence ratio higher than 1.5. Although such 
operating regimes are outside the usual HCCI operating conditions, 
this finding indicated that at low temperature-high equivalence 
ratio combustion conditions HC emissions can be reduced by sub- 
stituting DME for diesel fuel. The optimal operating range for the 
best ISFC or imep is at low peak cycle temperature and below stoi- 
chiometric for both fuels, because ignition is near TDC for that area. 
Total heat release of DME is close to that of n-heptane in this study. 
The optimal operating range, which was based on a merit func- 
tion including emissions and performance, was between 1700 and 
1850K peak cycle temperature and below 0.5 equivalence ratio for 
DME HCCI engines. 

Kong [142] studied numerically and experimentally the dual 
fuel DME-methane oxidation in an HCCI engine. The researcher 
used an improved version of the CFD code KIVA-3V [121]. The 
CHEMKIN chemistry solver [94] was integrated into the KIVA-3V 
code for solving the chemistry during multi-dimensional engine 
simulations. A detailed reaction mechanism for DME [53-55] was 
used to simulate fuel chemistry, consisting of 79 species and 351 
reactions. This mechanism was enriched by NO,-related reactions 
obtained from a NO mechanism, which was derived from the Gas 
Research Institute (GRI) NO mechanism [84]. The final scheme con- 
sisted of 83 species and 360 reactions. These reactions included the 
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oxidation chemistry of methane, which was used to simulate nat- 
ural gas. Turbulence was modeled via the RNG k-e model and a 
new wall function was constructed to calculate the wall heat flux, 
which accounted for the variation of the gas density and the turbu- 
lent Prandtl number in the boundary layer assuming quasi-steady 
conditions. 

Using this model, the researcher found that as DME concen- 
tration decreased ignition delayed and combustion duration was 
prolonged. In contrast to methane, DME experienced two-stage 
ignition and was depleted before any significant oxidation of 
methane. During the first stage of heat release, a noticeable amount 
of CO was formed, which was later increased during combustion 
and decreased rapidly only after both methane and DME were 
nearly consumed. The combustion of methane was induced by the 
high temperatures resulting from DME combustion. NO, emissions 
were much lower than those in conventional diesel engines, due 
to the low combustion temperature. It was found that there was a 
minimum requirement in DME concentrations for HCCI combus- 
tion regardless of natural gas concentrations, in order to obtain 
stable combustion. Moreover, the operating limits - defined by 
misfire and knocking combustion - narrowed as the methane con- 
centration increased, implying that high methane concentrations 
induce combustion instability under the conditions examined. This 
was probably the result of the rapid combustion of methane, which 
was triggered by the DME oxidation and prevailed as the amount 
of methane increased. 

Yu et al. [143] used large eddy simulation (LES) to simulate 
ethanol HCCI combustion, in which the ethanol was port injected. 
The main scope of the study was to investigate the role of turbu- 
lence in the HCCI combustion process. The chemical mechanism 
used was that of Marinov [60]. The LES simulations were conducted 
for two different piston configurations, i.e. a disc shaped combus- 
tion chamber and a bowl in piston configuration, while maintaining 
the compression ratio. A grid consisting of 128 x 128 x 128 mesh 
cells was used, with size about 1 mm in the radial direction and 
0.32-1.5 mm in the axial direction for the square bowl-in-piston 
case and 0.075-1.3 mm for the disc-shaped combustion chamber 
case. This mesh size was of the order of Taylor microscale and at 
TDC close to the Kolmogorov microscale. Eight parallel processors 
and 48 h were necessary to yield a full cycle simulation for a four- 
stroke engine. The simulation results, which were validated against 
experiments, showed that the piston geometry affected the HCCI 
combustion process; the square bowl-in-piston design had a much 
slower combustion. Although the levels of turbulence in terms of 
the velocity components were similar, the temperature inhomo- 
geneities were quite different between the two configurations. This 
was the reason for the longer combustion duration in the case of 
square bowl-in-piston design. 

A similar study was later conducted by Joelsson et al. [144]. In 
this study, the researchers compared numerically and experimen- 
tally the ethanol HCCI combustion in a metal and an optical engine. 
The engines were essentially identical, the main difference being 
the piston material; a metal piston was used in the metal engine 
and a quartz piston head in the optical engine. It was expected that 
the quartz piston engine would operate at a higher piston tempera- 
ture, since the quartz piston was air cooled, in contrast to the metal 
piston that was oil cooled. The scope of the study was to investi- 
gate the resulting thermal stratification differences and their effect 
on ethanol HCCI combustion. The same LES model was used in this 
study, as in the previous one. The expected difference in the piston 
temperatures was accounted for by increasing the piston surface 
temperature in the optical HCCI engine simulation. The simulation 
results showed that the two engines had minor differences in the 
induced turbulence field. However, the thermal stratification was 
much higher in the optical engine, due to the hotter quartz piston 
temperature and the lower intake temperature used in the optical 


engine. These factors prolonged the combustion duration, thereby 
decreasing the maximum pressure rise rate and the peak pressure. 

Viggiano and Magi [145] used a multidimensional CFD model 
coupled with chemical kinetics to study ethanol HCCI combustion. 
The model used was the REC-200, which solves the 3D Reynolds 
averaged Navier-Stokes equations for transient, two-phase, tur- 
bulent, chemically reactive flows with sprays [146]. Turbulence 
was modeled with a two-equation k-e model, and the heat and 
momentum fluxes were calculated using the wall function model of 
Launder and Spalding [147]. Combustion was modeled by coupling 
the turbulent flow to a detailed reaction mechanism. The approach 
used for the coupling was that of Kong and Reitz [148]. The influence 
of turbulence on the reaction rates was accounted for by multiply- 
ing the latter with a correction factor, which accounted for both 
the kinetic and turbulence timescales. The detailed reaction mech- 
anism used was that of Saxena and Williams [59], which included 
43 species and 235 reactions excluding NO, formation reactions. 
The structured grid of the 1-degree sector considered, consisted 
of 26 cells in the radial direction. In the axial direction, the cell 
number varied from 48 at BDC to 8 at TDC by removing cell rows 
during compression. The grid resolution near the wall was kept 
at 0.2 mm. After the validation of the CFD model against experi- 
ments was accomplished, the researchers investigated the effect 
of turbulent timescale, turbulent diffusivity, initial and wall tem- 
perature, and equivalence ratio on HCCI combustion. They found 
that increasing the turbulent timescale induced a reduction to the 
combustion rate and decreased the pressure rise rate. Increasing 
the initial diffusivity delayed ignition and deteriorated the com- 
bustion rate. Inlet and wall temperature affected combustion in a 
similar manner, with the effect of inlet temperature being more 
pronounced; specifically, increasing Tin or Ty advanced ignition, 
enhanced combustion, and reduced CO emissions. 

In Ref. [149] Liu and Karim examined numerically the effect of 
various operating parameters on hydrogen HCCI combustion. For 
the purposes of this investigation, the KIVA3 multi-dimension CFD 
code [121] was adopted with some modifications to make it suit- 
able for modeling the turbulent combustion processes. Detailed 
chemical kinetics of hydrogen oxidation was implemented into 
the code, including 20 reactions with additional reactions involv- 
ing nitrogen [71,72]. With this model, the researchers conducted 
a parametric study, the results of which showed that by increas- 
ing the swirl ratio heat transfer was enhanced, thereby reducing 
both the minimum and maximum temperature within the com- 
bustion chamber. As a result, ignition was retarded, combustion 
duration was prolonged, and the peak pressure and NO, emissions 
were reduced. An increase in engine speed increased the heat loss 
rate due to the induced increase of the swirl ratio, but the heat 
losses were reduced because the time available for heat transfer 
was decreased. The calculations showed that there was an opti- 
mum engine speed in terms of indicated thermal efficiency. Initial 
temperature and compression ratio affected significantly the onset 
of combustion as well as the combustion duration. 

Cotino and Jeanmart [150] studied numerically the running zone 
of ethyl acetate and compare it with the running zone of isooctane. 
Ethyl acetate was chosen, since it is a product of the acidogenic 
fermentation process used to produce biodiesel. This molecule was 
the basic compound and was used as a starting point to character- 
ize a group of esters produced. The model used in the numerical 
study was the commercial CFD FLUENT code and the computa- 
tional grid was obtained from the commercial software GAMBIT. 
The combustion mechanism used for ethyl acetate was a reduced 
one, consisting of 41 species and 203 reactions. This mechanism 
was derived from more detailed ones, namely the Gasnot et al. 
ethyl acetate submechanism [65], and the low carbon combustion 
mechanism of Miller and Melius [151]. Turbulence was modeled 
with the RNG k-& model and the molecular mixing at the subgrid 
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level via the eddy dissipation model (EDC). Except for swirl motion, 
3D effects were not included and an axisymmetric geometry was 
chosen, yielding 33,000 cells at IVC. 

The model was validated against experimental results using 
isooctane as the fuel, and was subsequently used with ethyl acetate 
to determine the ethyl acetate HCCI operating limits (running zone) 
at 1200 rpm. The upper load limit was determined by the maximum 
pressure rise rate, which was set at 5 bar/CAD. The lower limit cri- 
terion implied combustion instability and was determined taking 
into account the combustion timing, the EGR rate, the sensitivity of 
the combustion on the initial temperature, the compatibility of the 
final and initial temperatures, and the CO emissions. The results of 
the study showed that the ethyl acetate had a similar running zone 
range as isooctane, but it was shifted toward higher loads. More- 
over, ethyl acetate presented extended low temperature reactions 
compared to isooctane. 

Um and Park [152] analyzed the chemical effects of biodiesel 
on combustion and emissions of an HCCI engine fueled with 
diesel/biodiesel blends. For this purpose, they used a modified ver- 
sion ofthe KIVA-3V code [121] in a2D computational mesh, coupled 
with CHEMKIN chemistry solver [70]. The mechanism for biodiesel 
oxidation used was the one suggested by Brakora et al. [78], and 
was modified by Um and Park [153] to produce a reasonable igni- 
tion delay. The mechanism consisted of 156 reactions among 54 
species. Since it was found that neat methyl butanoate was not 
suitable as a biodiesel surrogate, due to its high oxygen content, it 
was assumed that 1 mol of biodiesel consisted of 1 mol of methyl 
butanoate and 2 mol of n-heptane. Neat heptane was also used as a 
diesel surrogate. A NO, mechanism, reduced from the Gas Research 
Institute (GRI) NO mechanism [84], was used for this study. It con- 
sisted of 12 reactions and 4 species (i.e., N, NO, NO2, and N20). 
Moreover, a soot formation and oxidation model was incorporated, 
since rich mixtures with equivalence ratio up to 2.0 were used. The 
high equivalence ratio, which is impractical in actual HCCI engines, 
is justified considering the nature of the study, which was to exam- 
ine the chemical properties of the biodiesel blends. Soot formation 
was calculated using a modified Hiroyasu model, which considers 
acetylene (C2H2) as a soot precursor [140], while soot oxidation 
was modeled using the Nagle-Strickland-Constable model [141 ]. 

The results of the study revealed that at stoichiometric equiva- 
lence ratio, increasing the biodiesel-to-diesel fuel ratio advanced 
ignition and decreased the total amount of heat released. This 
reduced the in-cylinder temperatures and induced a decrease in 
NO, emissions with an increase in CO and HC. At richer mixtures 
($» 1.0) increasing the biodiesel content increased CO, decreased 
NO, and HC emissions, and decreased Soot emissions compared to 
neat diesel fuel operation. Below stoichiometric conditions (4 < 1.0) 
CO and HC were formed in the cases where peak cycle temperatures 
were below 1400 K, NO, was produced at peak temperatures higher 
than 1800K, while Soot emissions were minimal. At rich mixtures 
significant amounts of CO and HC emissions were observed for most 
cases and soot emissions, which were higher for neat diesel, peak- 
ing in the cases with a peak cycle temperature of 1600 K. The NOx 
and CO emission maps, i.e. equivalence ratio versus peak tempera- 
ture maps, were similar for the diesel fuel and biodiesel cases. 


6. Conclusions 


Chemical kinetics are an essential part of biofueled HCCI com- 
bustion modeling, owing to the nature of HCCI combustion and 
the diversity of the biofuels’ molecular structure. The development 
of chemical reaction mechanisms is an active research area, from 
which biofuel HCCI research can benefit. Efforts should be focused 
on the consistency of the various mechanisms developed for use in 
HCCI engines. Special care should be exercised on the selection of 


a chemical mechanism, depending on the conditions under which 
its validation was conducted. These conditions are preferably 
high pressure and low equivalence ratio, since such conditions 
are encountered in HCCI combustion. It would also be advisable 
to implement more than one chemical mechanism for the same 
biofuel in any HCCI simulation model, and compare the results. 
This helps to avoid reaching conclusions (e.g. on autoignition, or 
emissions formation), which are the result of the specific chemical 
mechanism used rather than the HCCI combustion process. 

Single-zone models for biofuel HCCI combustion are very use- 
ful in identifying reaction paths and dominant chemical species 
during combustion. They have been used to provide insight into 
the behavior of various biofuels and blends thereof, and to validate 
the corresponding combustion mechanisms versus experiments or 
other more detailed mechanisms. Owing to the zero dimensional 
nature of the single-zone models, they have been used mainly to 
provide insight into the autoignition properties of the fuel and the 
formation of NO, emissions. Single-zone models also provide an 
indication on the effect of various configurations, such as nega- 
tive valve overlap, heavy EGR, fuel reforming, etc., on biofuel HCCI 
combustion. 

Multi-zone models with their added feature of accounting for 
mixture thermal and fuel stratification, provide better and more 
realistic estimates on combustion duration, peak combustion pres- 
sure and emissions formation. Multi-zone models for biofuel HCCI 
combustion have been used to study the effect of various biofuels 
and their blends on the in-cylinder thermal stratification, combus- 
tion duration, emissions formation, and the overall performance. 
Not only can they provide an insight into the formation of CO and 
HC emissions, but also to the specific species formed and consumed 
in the low temperature regions ofthe combustion chamber, thereby 
indicating the species expected to be found in the exhaust of a HCCI 
engine. This is an important feature, since the composition of HC 
emissions can vary substantially in the case of biofuels, due to their 
non-conventional molecular structure. 

Multi-dimensional models provide an even higher spatial res- 
olution of the combustion chamber than the multi-zone ones. Gas 
flow during the gas exchange process or within the combustion 
chamber can be modeled in more detail, providing a more realis- 
tic spatial resolution of the species evolution during combustion. 
Moreover, mass flow between hot and cold regions is described in 
detail, providing an insight into the emissions formation process. 
Especially in the case of direct injection biofueled HCCI combus- 
tion, the physical properties of the fuel must be described in detail; 
properties such as viscosity, volatility, etc., are essential to describe 
the evolution of the fuel spray into the combustion chamber. The 
importance of these properties cannot be overemphasized, since 
the spray evolution will determine not only the fuel stratification 
within the combustion chamber but also the fuel impingement on 
the cylinderlinerorthe piston crown. The description ofthe physic- 
ochemical properties of biofuels is an active area of research, which 
is bound to be further advanced. 
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